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In this study we address the question of the thermal description of collisions of classical clusters in 
the framework of classical N-body dynamics. We compare the results of systems in central collisions 
and those of thermalised systems with the same sizes and with the same available energies. The 
comparisons are made on size distributions, on total multiplicities distributions and on IMF (N > 3) 
multiplicities distributions. These distributions are found identical for available energies per particle 
below the energy of the least bound particle of the total system. They are notably different for higher 
available energies. These results may put a new light on the standard explanation for fragment 
formation in N-body collisions. 
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During the last two decades, strong experimental and 
theoretical efforts have been put on the interpretation 
of fragment formation in nucleus-nucleus collisions. This 
question is of great importance for the determination of 
the properties of the infinite nuclear matter. Experimen- 
tally a strong enhancement of fragment multiplicity with 
a charge higher or equal to 3 has been observed for ex- 
citation energies per nucleon greater than to roughly 3-4 
MeV. The usual way to theoretically describe such an 
enhancement is to assume a two step process: a dy- 
namical stage during which light charged particles are 
promptly emitted and hot sources are formed; and a de- 
cay stage during which the hot equilibrated sources de- 
cay through thermal emissions of particles and fragments. 
Many codes have been developed to describe these two 
phases, but the question of whether the fragments are 
produced during the first dynamical stage or during the 
thermal decay stage is not yet clearly fixed. Even if frag- 
ment size distributions have been successfully described 
by statistical fragmentation models, many experimental 
results indicate that the angular distributions and the 
kinematical properties of the fragments are strongly in- 
fluenced by the entrance channel. 

In this article, we will address the following question 
in the framework of the Classical N-Body Dynamics code 
(labeled CNBD): can the fragment production be de- 
scribed by the decay of a thermalised system? We hence 
have compared the results issued from the collisions of 
classical systems to the results issued from the decay of 
thermalised systems. In this article, we have limited our 
comparison to the most central collisions. 

First of all we briefly describe the CNBD code and the 
building of thermalised systems with it. Then, we will 
compare the total multiplicity distributions, the fragment 
multiplicity distributions and the fragment sizes distribu- 
tions. Finally we will draw conclusions. 

The dynamical code CNBD is described in details in 
references 0, Q . The evolution of each particle is driven 
by the Newtonian equations of motion. The two-body 
potential used to describe the interaction between the 
particles is a third degree polynomial. It has the basic 



properties of the Lennard- Jones potential 0,0, i-e. it 
is repulsive at short distance, attractive at medium dis- 
tance, and null at long distance. Since one wants to study 
the simplest case, neither long range repulsive interaction 
nor quantum corrections like a Pauli potential have been 
introduced @. A system is simply defined as a set of N 
particles. When the configuration of these particles min- 
imizes the total energy, the system is considered to be in 
its "ground state" . The dynamical evolutions of the par- 
ticles of colliding systems and the particles of thermalised 
systems are described with this code. 

In order to avoid any confusion with nuclear physics, 
the units used here are arbitrary. The distance will then 
be in Distance Simulation Units (D.S.U.), the energies in 
Energy Simulation Units (E.S.U.), the velocities in Ve- 
locity Simulation Units (D.S.U./T.S.U.) and the reaction 
time in Time Simulation Unit (T.S.U.). 

The "ground states" of systems made with CNBD are 
qualitatively close to that of nuclei Q ; and results of col- 
lisions between such classical clusters have already been 
shown in references 0, 0] . Central collisions are defined 
as those that have an impact parameter b lower than one 
tenth of the maximum impact parameter b max , which is 
the sum of the radii of the projectile and the target and 
the range of the two-body potential ( . < 0.1). For 
these collisions, the fragment (cluster with at least four 
particles) multiplicity has been found to strongly increase 
when the available energy per particle is above the en- 
ergy of the least bound particle E^ east Bound of the fused 
system. This energy has also been found to be the above 
limit of the thermal energy per particle that a free clus- 
ter can sustain . To perform the comparison of ther- 
malised systems with those in collision, we have build 
thermalised systems in the following way: i) the stable 
system of interest is put in a sphere with a radius R con f- 
Its center of mass coincides with the center of the sphere. 
Particles cannot escape this sphere because of a recall 
force derived from a quadratic potential, ii) The system 
receives a given amount of kinetic energy, iii) Then the 
particles move in the sphere according to the classical 
equation of motion until the system is thermalised. The 
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thermalisation is estimated by looking at the ratio < ^, c > , 
where tr c and < E c > are respectively the dispersion and 
the mean value of the kinetic energies of the particles. 
When this ratio becomes and then stays close to yj 3^-1 

with 10 -2 near, the system is considered as thermalised. 
This latter value is the expected value of the former ra- 
tio for an isolated system of size N for a classical system 
in the microcanonical ensemble |lCj. This procedure is 
equivalent to the procedures used in other works (see for 
example jHH). 

At this stage we could compare directly the results 
issued from the thermalised systems to those of the sys- 
tems in collision at the freeze-out time which is the time 
at which the fragments are well separated in the con- 
figuration space and do not interact anymore with each 
other. Unfortunately this time is extremely hard to de- 
termine for each collision. It has been more easier to let 
the thermalised freely decay during a time correspond- 
ing approximately to the difference between the ending 
time of the collision and a very rough determination of 
the freeze-out time. Hence, after the thermalisation time, 
the confining sphere is removed and the particles of the 
thermalised systems can move freely during roughly 170 
T.S.U.. This time does not need to be precisely defined 
since at large collision time, the emission rate is very low 
and the size distribution and the multiplicity distribution 
evolve very slowly. 

All along this process, the evolution of the classical 
particles in the system are followed by using the same 
algorithm as the one which is used for describing the sys- 
tems in collision, and the same two body interaction is 
used. This allow us to keep everything under control and 
to be sure that the similarities or the differences between 
these two systems will be only due to the different con- 
ditions of preparation of the systems. 

The comparisons have been performed for systems with 
26, 68, 100 and 200 particles. The results of a ther- 
malised system with N particles have been compared to 
the most central collisions ( , < 0.1) of a projectile 
with N/2 particles striking a target with N/2 particles. 
The same total energy in used in the two cases: the ex- 
citation energy of the fused system in the colliding case 
has been given to the thermalised systems. The compar- 
isons have been done for four excitation energies values: 
i) an excitation energy per particle E * /N below the 
binding energy of the least bound particle to the system 
ELeastBound ; h) an excitation energy per particle close to 
EheastBound] hi) an excitation energy per particle close to 
the binding energy per particle of the system Esind/N] 
iv) and an excitation energy per particle greater than 
Esind/N. For the thermalised systems, the influence of 
the size of the confining sphere has also been checked. 
Three effective densities have been studied: p = po 

{Rconf — Rmax(N)), p = po / 3 (R C onf = ^RtnaxiN)) 

and p = po/8 (R con f = 2R max (N)); R ma x{N) is the dis- 
tance between the centre of mass of the considered system 
and the particle which is the most distant from it. Since 
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FIG. 1: Normalized distributions of the total fragment mul- 
tiplicity for the system N = 68. The available energy in the 
center of mass increases from left to right and from top to 
bottom (see text). The full lines correspond to the central 
collisions of a projectile with 34 particles hitting a target with 
34 particles. The dashed lines correspond to the thermalised 
systems at normal density (p = po)- 



for the central collisions one hundred events are selected, 
one hundred thermalised events have been generated for 
each set of parameters (E * /N,N,p). 

The distributions of the total multiplicity, of the inter- 
mediate mass fragments multiplicity, and of the sizes of 
the final clusters are shown on figure ^ [21 an( i EH respec- 
tively. The intermediate mass fragments (labeled IMF) 
are defined as the fragments with at least four particles 
(N > 3). Theses distributions are shown for the system 
with N = 68 particles and at the four excitation energies. 
The dashed lines correspond to the thermalised case and 
the full lines to the central collision case. On each figure, 
the distributions corresponding to the low excitation en- 
ergy (E*/N < E LeastBo und) are displayed on the top left 
panel, those corresponding to E* /N close to EL eastBoun d 
are displayed on the top right panel, the distributions cor- 
responding to E * /N close to Esind/N are displayed on 
the bottom left panel and those corresponding to the high 
excitation energies (E* /N > Esind/N) are displayed on 
the bottom right panel. 

On figure ^ one can clearly see that at the low exci- 
tation energy (top left panel), the total fragment multi- 
plicity distributions for the central collisions and for the 
thermalised case are almost identical. For the highest 
energies, these distributions differ slightly at E * /N w 
ELeastBound (top right panel) and at E * /N ps E B ind/N 
(bottom left panel), the average multiplicity value ob- 
tained for the thermalised systems is greater than the 
average multiplicity obtained for central collisions. At 
the highest excitation energy (bottom right panel), the 
two distributions are clearly different: the thermalised 
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N 


E*/N 


Total Multiplicity 


IMF Multiplicity 


N of the biggest cluster 


Least Bound] 


(E.S.U.) 


C 


T 


C 


T 


C 


T 




37 


1.7 ± 0.6 


1.6 ± 0.6 


1 ± 


1 ± 


25.3 ± 0.6 


25.4 ± 0.6 


26 


60 


7.2 ± 1.0 


6.8 ± 1.3 


1.4 ± 0.5 


1 ± 


17.0 ± 2.9 


19.6 ± 1.4 


(50 E.S.U.) 


87 


12.0 ± 1.3 


12.4 ± 1.4 


2.1 ± 0.8 


1.5 ± 0.6 


7.5 ± 2.0 


10.1 ± 2.9 




107 


14.7 ± 1.9 


15.4 ± 1.7 


1.2 ± 0.7 


1.2 ± 0.6 


4.9 ± 1.4 


6.1 ± 2.1 




53 


6.2 ± 1.5 


5.9 ± 1.5 


1 ± 0.1 


1 ± 


62.5 ± 1.6 


62.9 ± 1.6 


68 


73 


15.4 ± 2.3 


16.8 ± 1.9 


1.7 ± 0.8 


1.1 ± 0.3 


43.7 ± 9.3 


50.2 ± 2.3 


(60 E.S.U.) 


108 


29.6 ± 2.4 


32.6 ± 2.6 


4.3 ± 1.2 


2.3 ± 1.1 


13.8 ± 3.1 


23.3 ± 5.5 




133 


36.3 ± 2.9 


40.4 ± 2.4 


3.9 ± 1.1 


2.9 ± 1.0 


8.7 ± 2.3 


10.7 ± 3.9 




42 


2.7± 1.1 


2.5 ± 1.1 


1 ± 


1 ± 


98.3 ± 1.2 


98.5 ± 1.2 


100 


72 


18.6 ± 2.9 


20.5 ± 2.2 


1.5 ± 0.7 


1.0 ± 0.2 


73.1 ± 11.9 


78.5 ± 2.6 


(65 E.S.U.) 


102 


35.1 ± 3.4 


41.3 ± 3.0 


5.0 ± 1.0 


2.4 ± 1.1 


23.4 ± 4.4 


43.4 ± 7.5 




132 


48.1 ± 3.6 


54.2 ± 3.2 


6.1 ± 1.2 


4.2 ± 1.5 


13.0 ± 3.6 


18.2 ± 6.1 




42 


2.3 ± 1.1 


2.3 ± 1.1 


1 ± 


1 ± 


198.7 ± 1.1 


198.7 ± 1.1 


200 


72 


24.6 ± 3.4 


27.6 ± 3.0 


1.1 ± 0.3 


1.1 ± 0.3 


174.1 ± 4.0 


171.1 ± 3.7 


(65 E.S.U.) 


102 


54.2 ± 4.6 


66.3 ± 4.1 


5.8 ± 1.1 


2.3 ± 1.0 


49.9 ± 9.2 


116.9 ± 8.6 




132 


81.5 ± 4.8 


94.7 ± 3.7 


10.2 ± 1.9 


6.2 ± 2.1 


24.5 ± 4.7 


46.4 ± 14.6 



TABLE I: Mean values of the total multiplicity ,of the IMF (iV > 3) multiplicity, and of the size of the biggest cluster for both 
central collisions (C) and thermalised systems (T) for different system sizes. The standard deviations (a) of these observables 
are indicated after the ± symbol. The density of the thermalised systems is equal to po. 
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FIG. 2: Normalized distributions of the IMF multiplicity 
(TV > 3) for the system N = 68. The available energy in 
the center of mass increases from left to right and from top 
to bottom (see text). The full lines correspond to the central 
collisions of a projectile with 34 particles hitting a target with 
34 particles. The dashed lines correspond to the thermalised 
systems at normal density (p = po). 



FIG. 3: Normalized distributions of the sizes of the final clus- 
ters for the system N = 68. The available energy in the center 
of mass increases from left to right and from top to bottom 
(see text). The full lines correspond to the central collisions 
of a projectile with 34 particles hitting a target with 34 parti- 
cles. The dashed lines correspond to the thermalised systems 
at normal density (p = po)- 



systems emit much more fragments than the central col- 
lisions. In the same way, the IMF multiplicity distribu- 
tions (figure |2J) are only similar at low excitation energy 
for the central collisions and the thermalised systems. 
For the two intermediate energies (E* /N « ELeastBound 
and E * /N EBind/N) the thermalised systems emit 
less IMF than the central collisions. This difference is 
slightly reduced at the highest excitation energy but the 



thermalised systems still produce less IMFs than the cen- 
tral collisions. 

As for the total multiplicity distributions and the IMF 
multiplicity distributions, the cluster size distributions 
(figure 01 for the thermalised systems and for the 
central collisions are almost identical for the lowest 
excitation energy. But they are very different at the 
two intermediate energies (E * /N « E^ east Bound and 
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E * /N w Esind/N): while for the central collisions the 
intermediate mass fragment area is well populated and 
no "fusion residue" can be clearly seen , the IMF area is 
empty at E * /N ps EheastBound and is barely populated 
at E * /N Esind/N for the thermalised systems. 
A clear contribution of a "fusion residue" is also seen 
at these two energies for the thermalised systems. At 
the highest excitation energy, the two size distribu- 
tions seem to be closer to each other, but the slopes 
of the distributions are different and the sizes of the 
biggest fragments are higher for the thermalised systems. 

From these distributions one can define two energy do- 
mains: one corresponding to E * /N values lower than 
EheastBound where the results from the thermalised sys- 
tems and from the colliding systems are almost identical, 
and one corresponding to E * /N values greater than 
EheastBound for which the distributions resulting for the 
two cases are clearly different from each other: while the 
thermalised systems emit mainly small fragments and 
keep a "residue", the central collisions produce mainly 
IMFs and no more "residue" is seen. This observation has 
to be linked to the conclusion made in the reference 0: 
the energy per particle that can be stored in a free cluster 
can not exceed the energy of its least bound particle. It is 
hence not surprising that the thermalised scenario does 
not correspond anymore to the central collisions when 
the excitation energy per nucleon of the fused system is 
greater than EheastBound- Another conclusion is that the 
thermalised systems can not produce a large amount of 
IMFs. To form an IMF, the particles have to be close in 
phase space. Such correlations are of course less prob- 
able in a completely thermalised case where the whole 
phase space is uniformly covered than in the dynamical 
case where only a small part of the available phase space 
is populated, due to the entrance channel effects. 

These observations seem to be independent of the sys- 
tem size. The table [I] summarizes the results obtained 
for the other systems sizes. The lines correspond to dif- 
ferent system sizes. In each cell, the average value of 
the observable of the corresponding column is indicated 
and is followed by its standard deviation. On the top 
of the column, the C letter corresponds to the central 
collisions and the T letter to the thermalised systems. 
The two energy domains can be clearly seen: the low en- 
ergy domain (E * /N < E Least Bound) where the results 
from the thermalised systems and the results from the 
central collisions coincide, and the high energy domain 
(E * /N > EL eas tBound) for which the results from the 
thermalised systems and the results from the central col- 
lisions differ significantly. 

This study shows that in CNBD the description of the 
central collisions as thermalised systems depends on the 
excitation energy per particle of the fused system. If this 
energy is lower than the binding energy of the least bound 
particle (EheastBound) of the fused system, the sizes dis- 



tributions, the multiplicity distributions and the IMF 
multiplicity distributions issued from central collisions 
are identical to those of the thermalised systems with the 
same sizes and with the same excitation energies per par- 
ticle. The thermalised systems decay through an evap- 
orative process which produces many small clusters and 
a big "evaporation residue" . This is consistent with the 
observation of a fusion/evaporation process observed at 
low energies in CNBD [j| . This behavior is qualitatively 
very similar to the behavior of Gemini code |3j which has 
been widely used to describe the thermal decay of hot nu- 
clei. When the excitation energy per particle is greater 
than E L eastBound-, the size distributions, the multiplicity 
distributions and the IMF multiplicity distributions of 
systems in central collisions are different from those ob- 
tained for the corresponding thermalised systems. While 
thermalised systems still produce mainly small clusters 
and an "evaporation residue", a large amount of IMFs 
are produced in central collisions and no more "residue" 
is seen. 

It is worth to notice that EheastBound seems to be a 
"key" energy for these classical systems as in references 
PI l^] : it is the fragmentation threshold for colliding clas- 
sical systems and the upper limit for the energy per parti- 
cle that can be stored in free clusters. This predominance 
of EheastBound, which can also be considered as a surface 
energy, shows that finite size effects play a major role 
for these classical clusters and can not be ignored for a 
full understanding of the reaction mechanisms of classical 
systems. 

Before extending these conclusions to nucleus-nucleus 
collisions, one has first to check the influence of a 
Coulomb-like long range repulsive interaction on these re- 
sults. The main deviation to these results may also come 
from the quantum mechanics. This last study is nowa- 
days unfortunately out of reach of such numerical simula- 
tions. One has also to remind that only central collisions 
have been studied in this paper and that for peripheral 
reactions possible geometrical effects could alter the ob- 
servations made here. This will be investigated in a forth- 
coming article. Nevertheless, providing that Coulomb- 
like effects and quantum mechanics do not qualitatively 
alter the conclusions made for these classical systems, 
this could shed a new light on the origins of the multi- 
fragmentation process in nucleus-nucleus collisions. Ad- 
ditionally, the exact meaning of the thermal description 
of the fragment production at high excitation energies 
(typically greater than 4-5 A.MeV) should be reconsid- 
ered. This could also lead to the re-interpretation of the 
liquid-gas phase transition signals seen in nucleus- nucleus 
collisions, like the negative heat capacity |IMIl9ll20| . the 
bimodality [II [M El [II and the Fisher's scaling [H . 
The occurence of such signals in the Classical N-Body 
Dynamics framework will be studied in forthcoming pa- 
pers. 
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